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I. INTRODUCTION 

According to Landau's Fermi liquid theorji^ low en- 
ergy properties of interacting fermion systems are simi- 
lar to those of an ideal Fermi gas. For many purposes 
interaction effects can be incorporated into renormalized 
parameters like the effective mass. Using the Landau 
theory one might expect that thermodynamic quantities 
like C(T)/T and xiT), where C(T) is the specific heat 
and x(^) is the spin susceptibility, have an expansion in 
powers of T^, as known from the standard Sommerfeld 
expansion for the ideal Fermi gas. 

This would mean extending the analogy too far, how- 
ever. In fact, it is known that interaction between 
fermions can induce so-called non-analytic corrections, 
that are absent for the ideal Fermi gas. The form of these 
corrections strongly depends on the dimensionality of the 
system. As a result of theoretical studies the leading tem- 
perature dependence of C(T)/T was found to be InT 
i n d = 3 and T in d = 2 at low temperatures, see Refs. 
HlUlStQ 10 and Refs. ITllllll3>.14.1,^ respectively. For 
the spin susceptibility a linear in T depende nce in d = 2 
was obtained in Refs. ll3ll4ll5ll6ll7ll8ll9ll2Clll2ll while 
for d = 3 a non-analytic In Q dependence of the 
wave- vector dependent spin susceptibilitjii^iiSiiLiSiS^ was 
found not to be paralleled by a similar temperature de- 
pendence. 

The model of a weakly interacting Fermi gas allows for 
a controlled perturbative expansion in the strength of the 
interaction potential. In such an approach the corrections 
cited above appear in the second order of perturbation 
theory. However, calculations in higher orders are diffi- 
cult and a full analysis of high order corrections has not 
been performed previously. 

Recently, a new low energy supersymmetric field the- 
ory for weakly interacting electrons was introduced in 
Refi. This is some kind of a bosonization that includes 
not only charge but also spin degrees of freedom. The 
scattering processes responsible for the non-analytic cor- 
rections are quasi one-dimensional in character, which 



leads to logarithmic contributions originating from inter- 
actions between spin excitations. As concerns the charge 
excitations in d > 1 studied previously within other 
bosonization schemea'^.i°i^3,^4i'^5,26,27,28,29,30,3i^ ^^^^^ con- 
tribution is less singular and does not lead to the effects 
discovered in Refi. The low energy theory of Refi re- 
sembles the supersymmetry approacb^SiS, well known in 
the theory of disordered systems and random matrix the- 
ory. 

In the model of interacting spin excitations, logarith- 
mic corrections to interaction amplitudes were found 
in any dimension, motivating a renormalization group 
study. The origin of the logarithms can be easily under- 
stood. Consider the interaction of two spin modes with 
propagators (vpyvp — iuj)~^, but opposite directions n, 
— n on the Fermi surface. Then one finds for perturba- 
tive corrections to the interaction amplitudes integrals of 
the type (compare Fig. ^ below) 



TY^Jdp 



1 



oc ln(eoo/T) 



(1.1) 



Eoo is the largest energy in the model. 

The integral over momenta transverse to n has to be 
regularized. Of course, in dimensions d > 1, the region 
of phase space for which the excitations move in almost 
antiparallel directions ni ^ — n2 is rather restricted and 
in general the logarithm is cut by max[r/eoo, |ni -I- n2|]. 
The question arises how the logarithms found at this level 
affect the temperature dependence of observable quanti- 
ties, since eventually one needs to perform a (weighted) 
average over all directions. To this end the specific heat 
was studied in Ref. It was found that C{T)/T gener- 
ally depends on an effective amplitude of backward scat- 
tering, that displays a complicated dependence on InT, 
in such a way that the results reduced to well known 
SCd=2 oc and SCd=3 oc T-^lnT when replacing the 
effective backscattering constants by the bare ones. 

The evident question arises whether similar corrections 
exist for the spin susceptibility. We will show in this 
paper that this is really so. To this end we generalize 
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the formahsm introduced in RefjlJ to include the exter- 
nal magnetic field. A convenient diagrammatic represen- 
tation makes apparent the relation to conventional dia- 
grammatic approaches. It turns out that, unlike for the 
thermodynamic potential, the contributions to the spin 
susceptibility in c? > 1 are determined by all the renor- 
malized interaction amplitudes of the model in various 
combinations. In d = 2 we find non-analytic contribu- 
tions Sx{T) oc T^l{T) and 7b (T) is an effective backscat- 
tering constant that includes logarithmic corrections. In 
d = 3 we confirm the absence of In T terms in the sec- 
ond order in the interaction. This does however not mean 
the absence of non-analytic corrections in three dimen- 
sions. Here we sum leading logarithmic corrections to the 
behavior and come to the result 5x{T) oc T^^^{T). 
Within the same formalism the one-dimensional case can 
be also considered and we reproduce the temperature de- 
pendence first obtained by Dzyaloshinskii and LarkinS^. 
Although the one-dimensional case is not the main focus 
of the approach, we consider the result obtained as an 
important check of the overall consistency. 

The paper is organized as follows. In Sec. we in- 
troduce the model that serves as a starting point for our 
subsequent analysis. In Sec. IIIII we decouple the inter- 
action part via a Hubbard-Stratonovich transformation 
in both the charge and spin channel and reformulate the 
partition function in terms of charge and spin excita- 
tions. Charge and the spin excitations decouple from 
each other and calculating the spin susceptibility we can 
concentrate on the spin sector. A representation for the 
partition function in this sector is derived in Sec. lIVI us- 
ing the supersymmetry technique. A low energy effective 
action is obtained and rules of the perturbation theory 
are introduced as well as a convenient diagrammatic rep- 
resentation. In Sec. 0we analyze corrections to the mag- 
netic field vertices in a renormalization group scheme. In 
Sec. IVII we study in our formalism temperature depen- 
dence of the spin susceptibility in one spatial dimension. 
Then we turn to the calculation of non-analytic correc- 
tions to the temperature dependent spin susceptibility in 
two and three dimensions in Sec. IVIII before concluding 
with a discussion of our results in Sec. IVIIII 



II. MODEL 

We introduce our model by specifying the partition 
function in the imaginary time formalisr" ^^'^^ 



Z = J V{x\x) exp(-5). 



(2.1) 



The fermionic fields x, x* depend on coordinates r and 
imaginary time r and carry a spin label a. They obey the 
antiperiodic boundary conditions Xo-lr, t) = — x^lriT + 
/3), where [3 = 1/T and T is the temperature. 
We write the action S as the sum of three parts 



5o describes free motion, Sb stands for the coupling of 
the spin to an external field b and Sint is the interaction 
of fermions, 



5o 



dxxli^) -~9r-HoXa{x) 



Sb = J dx x*a{x)bcr„„,Xa'ix), 



(2.3) 
(2.4) 



(2.5) 

Here and in the following summation over repeated spin 
indices is implied and v{x — x') = V{r — r')S{T ~ r'), 
where V{r — r') is the interaction potential. We use the 
notation 



(r,r) 



da; = 



dr / dr 



(2.6) 



and Hq = ^ — /J,, where p is the momentum operator 
and /i is the chemical potential. The standard relation 
between the partition function Z and the thermodynamic 
potential f2 is 



n = -TlnZ. 



(2.7) 



The spin susceptibility is then obtained as the second 
derivative with respect to b 



Xix-x') = - 



s'^n[b] 



5b{x)5b{x') 



(2.8) 



6=0 



We will mostly be interested in the static spin suscep- 
tibility Xs for a spatially homogeneous external field, in 
which case after Fourier transform the limit of vanishing 
frequencies should be taken before taking the external 
wave- vector to zero. 

Let us also introduce the following convention for the 
integration over momenta 



dq 



d'^q 
(2^ 



(2.9) 



S — Si) -\- Sb + Si, 



(2.2) 



in dimension d. For the sake of notational convenience 
we sometimes write J^, J^, J^, symbolizing J dq, J dr 

and dr respectively. / dh stands for the integration 
over the solid angle normalized to unity. 



III. DERIVATION OF THE BOSONIZED 
ACTION 

In this section we present the derivation of the model 
that will be used for the further analysis of the spin sus- 
ceptibility. It describes low lying charge and spin exci- 
tations in the system. A derivation in the absence of 
external sources has been presented in Ref. ^ Here we 
include the coupling to the magnetic field, so that we will 
mainly focus on the changes introduced by the magnetic 
field. 
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A. Decoupling into slow pairs 

For the interaction part we perform decoupling in two 
different channels by singling out slow pairs in the fol- 
lowing way 

Sint — *■ Sint — Sint,l + Sint,2 > (3-1) 

Sint,i = J dpidp2dqV{q} (3.2) 
Sint,2 = - J dpidp2dq V{pi - P2 - q) (3.3) 

^X*a{Pl)X<7'{Pl - q)X*a'iP2)Xa{P2 + q) ■ 

Here we denoted 



J dp.^TY, J rfp. 



(3.4) 



where the sum goes over fermionic frequencies e„ 
7rT(2ni + 1) and 



(3.5) 



where the sum goes over bosonic frequencies f2„ = 27rnT, 
further p, = (pi,e„J and q = (q, f7„). 

The cutoff function /, introduced in Eq. 1)3. 5|l . is de- 
fined as follows 



/(P) ^ foipro), P = IpI 



(3.6) 



where /o(i) has the properties /o(0) = 1 and f{t) 
smoothly for t oo. The function / has been introduced 
to avoid double-counting when singling out regions of 
small momentum transfer in Eqs. H3.1|l - H3.3|l . since with- 
out the cutoff both Sint,i and Sint,2 would identically 
reproduce the original Smt- Accordingly, A:c = r,^^ is a 
momentum cutoff that is much smaller than the Fermi 
momentum kc ^ pf, but much larger than typical mo- 
menta for the excitations of the low-energy theory that 
we wish to construct. We denote also the cutoff energy 
Eoo = ^'o'^^-F- Additional decoupling in the Cooper chan- 
nel is not included, since this would amount to overcount- 
ing of relevant scattering processes (compare the related 
discussion in Ref. Ill, Sec. II A). 

For a short range potential we can further simplify our 
considerations by setting V2 = ^(q Pf)- Since impor- 
tant momenta are close to the Fermi surface we can write 
14 (6*12) = l^(pi-P2) = y(2posin(^)), where 6*12 is the 
angle between momenta pi and P2, Ou — PiP2- 

For the further development of the theory it will be 
crucial to separate explicitly interactions in the triplet 
and singlet channel. 
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1, 



The action separates into a charge and a spin sector, 

Sint = Sint,s + Sint,t , (3-8) 

Sint.s = ^ J dpidp2dq p{pi,-q)Vs{9i2)p{p2,q) , 

Sint,t = J dpidp2dqS{pi,-q)Vt{9i2)S{p2,q) : 
where the charge p{p, q) and spin densities S{p, q) are 

pip,q) - ^^{p-l)^{p + l) ^ (3.9) 

S(P,9) = Xt(p-|)^x(p+|) , (3.10) 

and we turned to a spinor notation x = iXhXl)- 

Finally, one may decouple the interaction term Sint us- 
ing a Hubbard-Stratonovich transformation with a field 
4>nix) = i(pnix) + <Thn{x). Hcrc (pnix) and h.n{x) are 
real bosonic fields, so that 0n(r, r) = 0n(r,r -I- /?) and 
n is the direction of momentum p on the Fermi surface, 
n = p/|p|. The result is the following representation of 
the partition function 

2=-^ I Vcl)WMWt[h]Z[h,h,^] . (3.11) 

The weight functions Ws, Wt are shown below in 
Eqs. H3.18|l . (|3.19|) and A/" is a simple normalization con- 
stant that will not be displayed from now on. The parti- 
tion function Z[h, h, ip] describes the fermion motion for 
fixed configuration of fields b, h, ip, 

Z[h,h,p]^ J D{x\x) exp(-5e//[b,h,^]). (3.12) 

where the effective action Seff has the form 

S'e//[b, h, ip]^ So+ 5h[b] + j dpdridr2 



X'''(ri,r)0„ 



ri + 1'2 



x(r2,T)e^P('-i-'-^) .(3.13) 



By comparing Eq. (|2.4|) with Eq. (|3.13|) we observe the 
following simple relation 



5e//[b,h, (/?] =Seff[0,h + h,p], 



(3.14) 



which enables us to remove the field b from Sef / by a 
shift in h at the expense of changing the weight Wt ac- 
cordingly. 

Now we can write down a representation of the parti- 
tion function in the presence of the magnetic field as a 
weighted integral over field configurations 

Z = Jvcj) W,[p]Wt[h ~ h]Z[cj)] , (3.15) 



where 



14(^i2) = 1/2-itV^i(^i2), Vt{0i2) ^ -Vi{ei2) (3.7) 



m = J D{x*,x) exp(-5e//[0]) (3.16) 
5e//[0] = 5o + y"dpdridr2e^P("-^-"-^) 

Xt(ri,r)<^„ (''-l±Il\x{v2,r).{i.n) 
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The weights W^s[v3] and Wt[h] are 

Vl^s[(/?] — exp ~ ^ / dhidh2dqdT 



enters at coinciding points and this is why it is advanta- 
geous to Fourier transform with respect to the difference 
of coordinates. 



^*nA^,r) V-\0u,cD ^„.(q,T)] ,(3.18) G(x,xV) = / (dp) e^P^"^-""') Gp (^l±^,r,rV) -(3.24) 



- / dhidh2dqdT 



Wt[h] = exp 

h],,(q,r) K"'(^i2,q) h„, (q, r) ,(3.19) 



Eqs. (|3.1()I3.19|I represent the final result of this sub- 
section. 

The model is still quite general and in order to make 
progress further approximations have to be introduced, 
where the focus will be on calculation of Z[(j>]. Since it 
was possible to remove the field b from this term, the 
further derivation of the theory parallels that of Ref. Q 
up to the point, where the weighted integral over the 
field configurations of (j> is performed. In the following 
we will outline the main steps here in order to introduce 
our notations and prepare the subsequent discussion of 
the model. 



B. Bosonized action 



In Sec. IIIIB li the derivation of the representation for 
Z[(j)\ in the so-called quasi-classical approximation is out- 
lined following Ref. Q. Sec. IIIIB 21 deals with some addi- 
tional contribution from the high energy sector and the 
resulting change of weights. 



1. Charge and spin modes 

Integration over the fields XiX* the expression, 
Eq. ifTTEIl . for Z[(j)] results in 



= exp (Tr In (-9^-71:, 



(3.20) 



where the symbol Tr includes integration over coordi- 
nates as well as trace tr^ in spin space. The operator $ 
acts in accordance with Eq. H3.17|l . In the next step, we 
use a standard trick introducing an auxiliary integration 
over parameter u, which enables one to formally avoid ex- 
panding the logarithm while keeping track of appropriate 
combinatorial factors. Instead, the Green's function for 
fixed field configuration G (x, x'\u(f)) comes into play. 



exp 



Tr 



dudu In ( —dr — Hq + u$ 



-iTt / du 
'o 



(<Ig) {x,x\u(j)) 



(3.21) 
(3.22) 



m 

= exp 
The Green's function, 

G{x,x'\<p) = ^/ D{x,X*) X{x)x\x') e-5=/^ 1^1(3.23) 



Using the fact that Gp is sharply peaked at the Fermi 
surface one splits J dp ~ J dn J d^p, where ^p = 
p^/2m — /i and v is the single particle density of states 
at the Fermi surface per spin direction. After integration 
over ^p the quasiclassical Green's function 

5„(r,T,TV) = ^ J dCpGp(r,T,T'|0) (3.25) 

enters the expression 

Z[^]/Z[0] (3.26) 

= exp ^-iTTi^ J J du tT[(j)n{x)gn{r,T,T\u(j))]J . 

One of the main results of Ref i which we only cite here, 
is a set of decoupled differential equations for the charge p 
and spin S components of g^ix, r, r) in the decomposition 



namely 



where 



ingnix, r, r) = ipn{x) + S„(a:)cr, 

Ln.uSn{x,u) = -udrhnix) , 
Ln,OPn{x,u) = -udr(p{x) , 



-dr -\- ivpnX/ + 2iuhj, 



(3.27) 

(3.28) 
(3.29) 

(3.30) 



hn{x) is a matrix in the spin space with components hij = 
— eyfehfc, so that ft,S = h x S. 

This result was obtained with the help of a generalized 
Schwinger ansatz'^^ for g^ 

g„(r,T,T'|</)) =T„(r,T)5o(r-T')r„-i(r,r'), (3.31) 

where go is the Green's function for free fermions (0 = 0). 
It was further assumed that (f> varied smoothly on the 
scale of the Fermi wavelength Xp — pj,^. 

Finally, substituting the decomposition Eq. (|3.27|l into 
Eq. 1(3231, one finds 



Z[<j>] = Z[0]Zp[^]Zs[h] 



(3.32) 



where 



Zp[(p] = exp J dudxdh (pn{x)pn{x)J , (3.33) 
Zs[h] = exp ^— 21/ J dudxdh hn{x)Sn{x)^ ,(3.34) 



and S and p fulfill the differential equations Eqs. 

The equation for p is readily solved using a 
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Fourier transform but we concentrate in the folfowing 
on the spin sector here. 
Generally, one finds 



Z = Z[0]ZpZs 



(3.35) 



where 



Z„ 



J V^WMZM , (3.36) 
= j VhWt[\i-h]Zs[h\. (3.37) 
Our main task is to calculate the partition function Zs- 



2. Weight functions 

As is well known, the quasiclassical approximation 
does not capture contributions that originate from fre- 
quencies of the order of the Fermi energy sp- In the limit 
of weak interactions considered here we may incorporate 
such contributions into the model by replacing weights 
Ws, Wt of Eqs. (rnm - irnm by weights Ws, Wt 



Wt[h] 
where 



exp 
exp 



- / dxd-n. Lpri{x) 



dxdfi hn(a;) 



1 + 2iyVs 
and we adopted the convention 



{x,n) ;3.38) 
(x,n)b.39) 



1 - 2iyVt 



(3.40) 



fp {x,n)^ / dri/(r - ri)pn(ri,r) , (3.41) 



(3.42) 



Vt^sb (x,ni) = J dn2Vt^sidi2)bn2{x) 



The above argument is valid, however, only for the 
quadratic in h part of the weight function Wt entering 
Eq. (|3.37() . It docs not hold for the part containing b, 
since the Green's function, for which the quasiclassical 
approximation was used, depends only on the field h. 
Therefore, writing Eq. (|3.37|) we should use the following 
form for the weight Wt [h — b] 



Wt[h-b] ^ >Vt[h,b] = Wt[h] 
X exp 1 dx h{x) 



(x) 



xexp j dxdnh(x) Vf'^'h {x,n) .(3.43) 

Here and in what follows we omit the cutoff function /, 
whenever the momentum is determined by the external 



field b. This cannot change results since we are interested 
only in small external momenta, |g| ^ f^^i which allows 
us to put /(q) = 1. 

It is clear from Eq. H3.37|l that the magnetic field cou- 
ples only to the spin degrees of freedom and the charge 
sector does not play any role for the spin susceptibil- 
ity. Therefore, we can concentrate on the spin sector, 
described by Eq. H3.37|l with the weight Wt [h — b] deter- 
mined by Eq. H3.43|l . 



IV. SUPERSYMMETRIC REPRESENTATION 



In Sec. rrvn we derive a representation for Zg, 
Eq. H3.37|l . in terms of a functional integral over super- 
fields. A detailed derivation has been presented in RefQ. 
Therefore we only highlight the main ideas here and rel- 
egate more technical details of the construction of the 
model to Appendix ^ In Sec. IIV Bl we collect the rel- 
evant definitions of supervectors and supermatrices that 
enter the final model. This model is then presented in 
Sec. IIV CJl rules of the perturbation theory are formulated 
in Sec. IIV Dl and a convenient diagrammatic representa- 
tion is introduced in Sec. IIV El 



Zs as an integral over supervectors 



Using Eq. (|334jl and Eqs. (jX2S|) . (|XgU|) one arrives at 
the following form for the partition function Zg [h] 

Z,[h] = (4.1) 

exp I 21/ I dudxdn hn{x) uL^\j^drhn (x) 



If one could find L^\^ exactly for all u 0, the problem 
would be solved. However, since this is hardly possible 
for an arbitrary u and h, we have to resort to some ap- 
proximation scheme. For this purpose it is advantageous 
to reexpress L~\ in terms of a Gaussian functional inte- 
gral. Using either bosonic (complex) or fermionic (Grass- 
mann) fields separately one would have to deal with an 
h-dependent normalization factor of the Gaussian inte- 
gral, which is inconvenient. This complication can be 
avoided by introducing an integral that includes both 
bosonic and fermionic variables on equal footing, as it 
has been used for a long time in the theory of disordered 
systems, where the technique is known as the supersym- 
metry method^^-^^ . In the context of the present problem 
this approach has been introduced in Ref. 1. 

When using the Gaussian functional integration one 
should be careful, however, since the operator L is not 
hermitian. In particular, the sign of the h-dependent 
term is not fixed and thus the convergence of the Gaus- 
sian integral over bosonic variables, for which one re- 
quires a positive real part of the kernel, is not easily in- 
sured. Fortunately, it is known how to overcome this 
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difRculty'^^. One can construct from the operator L Her- 
mitian operators L' = {L + D)/2 and L" = -i{L-D)/2 
and arrange them into a new hermitian matrix operator 



M = 



L' iL" 

-iL" -L' 



(4.2) 



H 



The corresponding vector space wih be called "Hermi- 
tized" or H- space. One can reconstruct by sum- 
ming certain matrix elements of the inverse of as 
was shown in Ref. 

The implementation of the ideas presented above leads 
to the following identity 



Zs[h] = cxp (,y f T^{X)nx\,:Fi,{X'] 
\ Jxx' 

Here we use the collective variables 

X — [x,z), X = (r, t), z ~ {u, n) 
and the integration measure is specified as 



dX = I = 

X 



dxdz, 



dz = 



du / dn 



(4.3) 



(4.4) 



(4.5) 



where J dx has been introduced in Eq. H2.6|l . 

To make contact with the previous discussion, we note 
that (A7Y)~i corresponds to of Eq. where A 

is some constant matrix introduced below. J-"h and !F\i 
are linear in h and their role is merely to select relevant 
components of TL~^ , the sum of which gives . The 
operator Ti-^x' '^^^ be written in terms of a Gaussian 
functional integral as follows 



.—nx\>-{Mx') 



(4.6) 



= j V{'^, ^) VxV'x' exp (^2iiy J + iSA) ^x 



where ip and ip = ip'^A are supervectors (see below). They 
contain both complex and Grassmann variables on equal 
footing, which leads to the simple normalization of the 
Gaussian integral 



P(V', i^) exp \^2iiy J ■ipx{'H + iSA) tpx j ■ (4.7) 

The identities, Eqs. (|4.3() . (|4.6|l . are the basic build- 
ing blocks for the derivation of the model we want to 
work with. In the next section we will define all quanti- 
ties involved in more detail, making it possible to verify 
Eq. H4.3|l by explicit computation. Construction of the 
theory that follows Ref 1 and adopts the notations used 
in this paper is included in Appendix^for the interested 
reader. 

It may be worth making a comment concerning the 
angular integration J dn. It will turn out later that for 
our purposes the most important scattering processes are 
forward and backward scattering. By forward scattering 



we mean scattering processes, in which both the incom- 
ing and outgoing fermions have almost parallel momenta. 
Backward scattering refers to a process, in which both 
incoming and outgoing momenta for each fermion arc al- 
most anti-parallel to each other, while incoming and out- 
going momenta of different fermions are almost parallel. 

It is then convenient to split the angular integration 
/ dh into two half-spheres, one half-sphere contains "left- 
movers" the other one " right movers" . The arbitrariness 
involved in fixing the boundary in dimension d > 1 will 
not become important due to the quasi-one-dimensional 
character of the relevant scattering processes. When sep- 
arating sectors of left and right moving particles it is then 
only necessary to integrate over one half sphere and we 
denote this angular integration as J dn with normaliza- 
tion J dn — 1/2. As an example, the angular integration 
in Eq. (|4.1() is now written as 

J dnh^L-'drhn (4.8) 

= J dn h„L„^9^h„ + h_„LI^(9^h_„ . 

Let us perform two more manipulations to arrive at a 
form, where only the averaging with weight Wt remains 
to be done. Starting from Eq. 14.7(1 one may verify by 
shifting fields ip, 'tp, that 

Z,[h] = exp(^j^J^^^ 'K{X)Hx.x'Fh{X')^ (4.9) 

2?(V','0) exp ^2w J V'x(W + «^A)Vx 

X exp (^V^tiyJ^ (Fi,{X)i^x +i^xFh{X))^ . 

Just as L, 7i contains a part Tih, that is linear in the 
field h. We split off this part by writing 



7i = "Hq + Ti-h 



(4.10) 



The final form of the model is then obtained by averaging 
the h-dependent part of Zg [h] with the weight Wt [h — b] 
(compare with Eq. H3.37|) ). 

Zs = (4.11) 



V{^P, i,) exp ( 2iv 1^ (^0 + iSA) iPx ] B[^, ^j, b] 



B[i^, ^A, b] = / Ph >Vt[h - b] exp 



2iv / ipxT^h'tpx 

X 



X exp 




{Ti,{X)^x+^xW))^ (4.12) 



In the next subsection we will give the precise definition 
for V', ip, ^, ^ and H. Equations ((4.3|l and Eq. 1(4. 6|l are 
discussed in Appendix^ What remains then is to obtain 
the final model by explicitly computing B, Eq. ((4.12() . 
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B. Relevant supervectors and supermatrices 

1. Supervector ^ and its conjugation 

Let US first introduce the supervector if) depending on 
coordinates X = {x, z) (cf. Eq. (|4.4|l ). It has components 
in the sectors of left and right-moving particles labelled 
as n, the graded space of bosonic S and fermionic x vari- 
ables labelled as g, the Hermitized space labelled by H 
and the spin space labelled by s. An additional sector is 
introduced, which simplifies calculations with the model. 
It has been termed "electron-hole" eh space in Ref. Q 
and plays a similar role as the time-reversal sector in the 
cr-model description of disordered systemaSS, 



where 

V = 



si 

S2 



(4.14) 



H 



Both S' and are vectors in the spin space 



S ^ 



= I SI 

SI 



(4.15) 



The components Xx^XyjXl a-re anticomniuting (Grass- 
mann) fields. 

The conjugate vector ^ is defined as 



■0 = ip'^A, 



where 



A 



1 

-1 



(4.16) 



(4.17) 



H 



is the third Pauli matrix in the Hermitized space. An 
important symmetry that arises due to introducing the 
(eh) sector is 



where C is the following matrix 
and matrices Ci have structure in the eh sector. 



eh 

















C2 = 




J) 



(4.18) 



(4.19) 



(4.20) 



For Grassmann variables the convention (x*)* — —x is 
used. The conjugation of matrices is introduced as 



where the special transposition appropriate for 
supermatrices^ should be used. The important 
property 



ipAcj) = (j)Aip 



(4.22) 



where "0, </> are supervectors, is one of the main motiva- 
tions for introducing the (eh) sector. When calculating 
higher cumulants later using Wick's theorem the num- 
ber of contractions can be reduced considerably with the 
help of relation Eq. (|4.22|) . 



2. The matrix TL 

The matrix 7i is split into a h-dependent and a h- 
independent part 

H = Ha+Hh (4.23) 
Ho = -ivQTsY.sn'V - Aidr, Hh = -2zt3H„ 

Different constant matrices in this expression are 



and 



H„(a;) 



(4.25) 



hnix) 

h^nix) 
3. Vector 



Vector does not have the full symmetry in super- 
symmetric g-space. Instead, it projects onto the bosonic 
sector. The role of the fermionic fields in Eq. (|4.9() is 
only to provide the normalization. We present JFh as a 
product of an h-dependent and an h-independent part. 
The latter one is 



V 



1 
1 

-1 
1 



eh 



eh 



(4.26) 



The charge conjugated vector is JFq = (CJ-q) , where 
only C2 is effective when evaluating the right hand side. 
J-'h and J^h are then given as 



and Th — (CJ-h) , where 
dx{a) = 



1 

ii( aSr + (1 — a)iwonVE 




(4.27) 



(4.28) 



A = CA' C 



Tr<T 



(4.21) 



Here an additional parameter a has been introduced into 
the model and we will comment on it in the following 
subsection. 
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4- Parameter a and the weight Wt. 

In view of Eq. (|4.1() one would expect only to enter 
Eq. H4.28|l . corresponding to a = 1. Choosing different 
values of a, however, may be convenient as we will see 
below when studying the renormalization of the model. 
We will set a = 1/2 there, treating temporal and spa- 
tial derivatives in a symmetric way. For a ^ 1 rela- 
tion Eq. 14.4|l needs to be modified and this modification 
eventually changes the weight Wt- Any physical quan- 
tity calculated with the model is of course independent 
of the choice of a. As is shown in Appendix when in- 
troducing parameter a ^ 1, relation Eq. (|4.4|l takes the 
form 

Zs[h] = exp - a) J hlix)^ (4.29) 

X exp f -4^l/2 f F\^{X) {^x^x') FhA^') 
\ Jxx' / 

For a = 1 it coincides with Eq. 14.4|l and one can use 
the form for B given in Eq. (|4.12|l . For general a it 
seems natural to absorb the exponential in the first line 
of Eq. (|4.29|l , into the weight Wt and thus change B to 



2w / ^px'HhiJx 

X 



B[il), ^A, b] = j Vh Wt [h, b, a] exp 



X 



FhA^)^x+^xFi,AX) 



,(4.30) 



X exp 
where 

Wt[h,b,a] = Wt[h,b] e-"'i-")-fft.-^"("\ (4.31) 

Clearly, this change in Wt only affects the quadratic form 
in h but not the part containing b. Therefore, to make 
the change explicit we may write here 



>Vt[h,b = 0,a] 



(4.32) 



exp 



dx{dn)h^{x) tf^{a)h (x, n) 



where 



2ft(a) = / 



2iyVt 



1 - 2vaVf 



(4.33) 



The final step in the derivation of the model is the cal- 
culation of B in Eq. 102011 ■ 



C. Effective low energy theory 

From Eq. (|4.11|) together with B given in Eq. H4.3U|I we 
find 



Z, = j I?(^,V') exp (-^5, 



(4.34) 



Next we specify the different parts Si of the effective 
action. 

The interaction-independent part is 



So = -2w j dX ipx{T-(-o + iSA)ipx- 



(4.35) 



There are three different interaction vertices present in 
the theory 



52 = -2iiyY2x,j J dXdXi 



(4.36) 



{i^X,S'^3^jdxJ^o) ^X.Xi i^odxi^jTs^pXus) 

S3 = -4V^j/ Ay ssfjj J dXdXi (4.37) 



ii^xjUT^njijx^ls) T'x^x^ i^odx^TstPxui) 

54 = -4z/^ Ay £5/37e5/3i7i J dXdXi (4.38) 

{'^X,pUT3Uj^X,i) r^^xi (V'Xi,ft'^lT-3nj'(/'Xi,7i) 

Summation over spin indices is implied and we use the to- 
tally antisymmetric tensor EaP'y with £123 = 1. This part 
of the action would be sufficient for a calculation of the 
thermodynamic potential in the absence of a magnetic 
field and it coincides with the action written in RefQ. 
Here we used the notation 



rx,x'^l.{^') f{v-v')5{T-T') 



(4.39) 



and 



7i(fmi) = 

72(nnl) = 
Matrices II,- are 



vVt 



1 - 2vaVt ^ 

vVt 
1 - 2vaVt 



(iv5l)=70 (4.40) 
(nr^i)^70 (4.41) 



The form of 



Hi = 1, n2 - S3 



1 1 

1 -1 



(4.42) 



(4.43) 



was determined from the following identities 

= i (niAHi + n2An2) (4.44) 

+ = i (niAHi - n2An2) (4.45) 

Taking into account that relevant scattering events are 
quasi one-dimensional, as will be seen later, Eqs. 14.44|l . 

also explain the labelhng in Eqs. (|On|) . lUUTJ, 
where 71 is classified as forward scattering and 72 as 
backward scattering. 
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The presence of a magnetic field introduces three more 
terms, namely 

^60 = -yri j dx h'^{x) (4.46) 

5f,i = -2v./^ir, j dX hs{x) {^x^sradxTo) (4.47) 

Sb2 = 4i^ esi3j 7j J dX bs{x) {ipx^pUTsijx,^) (4.48) 

In these expressions 

V=-^-^, (4.49) 
1 — 2auVt 

where the bar in Vt means averaging over the full solid 
angle. We remind that rj, the interaction amplitudes Ti 
as well as dx, dx depend on parameter a introduced in 
14.281 Here we suppressed the label for the sake of brevity. 

D. Rules of Perturbation Theory 

1. Gaussian averages 

A perturbation theory can be set up in a standard 
way using a cumulant expansion and Wick's theorem. 
Gaussian averages are taken with respect to So and it 
is therefore convenient to work with the matrix Green's 
function 

g{Xi,X2) = -4i,y{^Px,(E>^x,)o (4.50) 
{{...)),= [ I?^(...)exp(-5o[V]) (4.51) 



Due to supersynimetry no normalization factor arises. G 
is a matrix in spin space but its spin structure is trivial 
and we denote 

Gal3iXi,X2) = Gniixi - 0:2 )<5a/3(5ni ,n2 (^("l " ""2). (4.52) 

The Fourier transform of G is introduced as 

e„i(a;i-X2)= (4.53) 
/ '^P ^ni(p,o;„)e^P^"-^-"-^)-^'^"(^^-"^). 

Here LUn — 2'iTTn are bosonic Matsubara frequencies and 

1 



t/„(p, w„) 



(4.54) 



Similar to the specific heat^, the temperature dependence 
of the susceptibility is determined by non-zero Matsub- 
ara frequencies only. Therefore, the term containing the 
infinitesimal i5 in the Green functions will not become im- 
portant in our calculations and will not be written from 
now on. The matrix Green function Gn is diagonal in 
spin space. 



In addition to averages of the type written in Eq. (|4.50|l 
one has to account for non-standard averages of the type 

('0Xi,aV'X2,/3)o and {il'x^.ai'x^,^)^ ^ 0' which arise 
due to the eh sector. It is, however, sufficient to work 
with Eq. H4.50|l . since expressions involving such non- 
standard averages can easily be transformed to a more 
standard form with the help of the relation tj^x-^ -^'4^X2 = 
xpX2^'*Pxi, which is valid for any supermatrix A. The 
charge conjugation operation A ^ A has been defined in 
Eq. ljOT|l . 

It is often convenient to work with a generalization of 
the trace operation tr, used for conventional matrices, to 
the so-called supertrace strt?^. It is defined as 



str 



a a 
P b 



— tra — tr& 



Two more useful relations are 



str(AB) = stT{BA) , 

^l,aM2,P = -Str(A (7/'2,aV'l,/3) 



(4.55) 



(4.56) 
(4.57) 



A, B are arbitrary supermatrices. It follows from the def- 
inition of the supertrace that for matrices Ai that have 
the full symmetry in ^-space, str(Ai) — 0. The ma- 
trix Green's function ^ is a particularly important ex- 
ample. One immediately concludes that the following re- 
lations (and straightforward generalizations thereof) hold 
for such matrices 

((V^AitA))^ = ((^Ait/.) (^A2iA))o = 0. (4.58) 

These important relations considerably reduce the num- 
ber of diagrams to be considered in the perturbation the- 
ory. 



E. Diagrammatic representation 

Figure n displays the building blocks that we will use 
for the diagrammatic representation of the perturbation 
theory. The three interaction vertices reflect the struc- 
ture of ^2, ^3 and 54. They are plotted in such a way 
that small momenta flow along the interaction line. The 
dotted lines symbolize the structure TqTq. These lines 
do carry neither momentum nor frequency. They are 
however convenient to make contact with conventional 
diagrams formulated in terms of electron Green's func- 
tions as we will discuss now. 

To this end it is instructive to make comparison with 
the original model. Since the charge channel has been 
separated, we need to consider fermions interacting in 
the triplet channel only, i.e. the free action Eq. (|2.3() and 
an interaction part of the form 



S --1 



S(pi,-9)rt(0i2)S(p2,9), (4.59) 



P1P29 
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FIG. 1: Basic building blocks of the perturbation theory. 

where we remind that S{p, q) — {p ^ f ) ""^ f) 
and we used four dimensional notation for momenta and 
energies. 

In contrast, the effective low energy model of Eq. (|4.34|) 
is formulated in terms of spin modes. The structure of the 
terms appearing after expansion of J hLhC^rh in Eq. (|4.1|) 
in the field h is the translation to the spin mode language 
of a closed Fermion loop with n > 2 fields h coupled to 
it. The expansion of Lh in powers of h is performed 
here assuming that Lh=o (wonp — zw)"^ describes free 
propagation of spin modes (compare to Sq of Eq. H4.35|l '). 
After integration over h, which reduces to contracting 
pairs of fields (hh), one obtains a theory of interacting 
spin modes. Not all fields h in / hLh^rh enter in an 
equivalent way, however, and this explains the presence 
of three different interaction terms in the model. When 
contracting two fields h that appear due to an expansion 
of Lh in h, one finds an interaction vertex of the type 
represented by ^4 (see Fig. ^ . If only one such field is 
involved one comes to iSs, otherwise to 

Let us summarize the discussion with the help of 
Figs, m O Here a closed loop of fermionic Green's func- 
tions (Fig.[2Jl as well as a particular diagram for the sus- 
ceptibility (Fig.|3J| are shown formulated first in terms of 
fermionic Green's functions and then also as a diagram 
for interacting spin modes in the effective theory. We see 
that one diagram of the conventional perturbation the- 
ory produces several diagrams of the expansion in the 
spin modes (we marked the corresponding interactions 
vertices by Si, i = 2, 3, 4). 

At first glance it looks as if our effective perturbation 
theory has become even more complicated than the origi- 
nal one. However, it is not so because its not sufficient to 
just write the diagrams. One should calculate them sin- 
gling out the most interesting low energy contributions. 
This singling out has already been performed when de- 



FIG. 2: A closed loop of fermionic Green's function contrasted 
with the spin mode representation. 




FIG. 3: A particular diagram for the susceptibility formulated 
in the fermionic and the spin mode representation. 



riving the effective theory for the spin excitations. A 
larger number of the diagrams in the new theory corre- 
sponds to different possibilities of obtaining low energy 
contributions when integrating in the diagrams of the 
conventional perturbation theory. 

The perturbative expansion obtained from the low 
energy effective action is equivalent to expressions ob- 
tained from Eqs. (|2.3() and (|4.59() after expanding ^p+q ~ 
^p + w_Fnq in the vicinity of the Fermi surface in each loop 
and subsequent integration in ^p. The present model, 
however, organizes the terms in a way, that is much more 
convenient for identifying the most important contribu- 
tions. 

Let us remark that the field b enters the diagrams in 
the same way as h and that for each loop there is one 
fixed angle n, which is as a direct consequence of the 
integration over ^p. 



F. Bare spin susceptibility 

The static spin susceptibility x at T = 0, 

X = 2vr^o.=u (4.60) 
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where rj is given by Eq. H4.49|l . can be obtained straight- 
forwardly from iSbo, Eq. (|4.46|l . for a = 1. For repre- 
sentations with an arbitrary a ^ I, the term Sbo alone 
does not provide the full answer. For general a one 
should consider the combination Sbo — 1/2 {{S'^i )) with 
Sbi from Eq. (|4.47|l . The average in this formula is to 
be taken with the full quadratic form Sq -1-52- This 
procedure effectively amounts to a ladder summation. 
In particular, we will later choose a = 1/2. In this 
case one arrives at Eq. (|4.6UI) with the help of the iden- 
tity 2i'T]i = 2vr]i/2 + 2vrj^^^{\ — i^Vtr]i/2)~^ , where the 
first term is obtained from Sbo and the second one from 
-l/2((5,^i)). 



V. RENORMALIZATION 

As sketched in the introduction, logarithmic correc- 
tions appear in the model and they can be summed in 
a renormalization scheme as it has been done for the 
interaction amplitudes in Ref. 1. The appearing loga- 
rithms and thus the renormalized amplitudes generally 
depend on the deviation in angles from the ideal forward 
or backward scattering case. It is therefore immediately 
clear that this renormalization scheme cannot include as 
effective charges physical quantities like the susceptibility 
that do not depend on such angles. This is why we con- 
sider the renormalization of interaction amplitudes and 
external vertices first and then include the renormalized 
values into a perturbation theory for the susceptibility. 



We decompose a supermatrix P in such a way that 



P = YA=i P" and 



[pi,n2] = o, {p\n3}-o (5.3) 

[p2,n2] = o, [p2,n3]-o (5.4) 

{p3,n2} = o, [p3,n3] = o (5.5) 

{p4,n2} = o {p^^3} = o (5.6) 



An explicit formula for P' can be given 



1 4 

P' ^ -Y,^^k^kP'Rk, K 



k=l 



1 1-1-1 

fc = I } _l 1 -1 1 (^-^^ 

1-1-1 1 



It is easily seen that 

4 

str[P,P,] - %str[p2], str[PQ] - ^ str[PkQk\- (5.. 



k=l 



The following useful relation can be checked by direct 
computation 



KikiK2k2 [^Lkik2BLkik2] 

ki ,k2 = l 

4 

4-5^,^,^ A,,fc [An^PEfc], (5.9) 



fc=i 

where Lfc^fe^ = Ilfc^IIfc^. 



A. Generalization 

During the process of the renormalization new terms 
may appear in the action. To consistently take such 
terms into account, they should be included into the 
model from the beginning. To this end a generalization 
of Eqs. i s intro duced in Sec. IV A 21 This 

step is prepared in Scc IV A 11 



1. Decomposition 

In Eqs. (|4.36|l - (|4.38(l the following decomposition of 
supermatrix P was used 

2 2 

p = ^p\ p' = i^A,fen,pnfe, (5.1) 

i=l k=l 

where [^^,£3] = and {yl^,S3} = for arbitrary ma- 
trices A. 

For a generalization we consider 4 supermatrices 11^ 
with [IlijIIj] = and Ilf = 1, namely 

Hi = l,n2 = S3,n3 = Air3,n4 = A1E3T3 (5.2) 



2. Generalized action 

After this preparation we introduce the generalized 
model. We start with ^4. 

Si = ~2iy ^ Ay esp^sspi-yi / dXdXi (5.10) 

The amplitudes Ti are 

f,(X,Xi) (5.11) 
= rj(nni; w, ui, (r - ri)^) /(r - ri) S{t - t') 

and by comparison with Eq. (|4.37(l one finds their bare 
values 

r.(^,.,ui;r-) = {;?{,^J (5.12) 

where = r — (nr) is the component of r transverse 
to n and one should keep in mind that important initial 
and final angles n, rii are almost parallel or almost an- 
tiparallel to each other. The initial values for F; do not 



12 



depend on u, ui and r-*- but develop such a dependence 
under renormalization. 

The generahzation of Eq. (|4.36|) reads 

S2=-iiyYl >^^J J dXdXi (5.13) 

ij cri,cT2— i 

In this formula J-± — t±Tq, T± = Tqt^ and t± — (1 ± 
T3)/2 are projection operators that change under charge 
conjugation as t± = Tip. In analogy to Eq. (|5.11|l we 
defined 

Af"^(X,Xi) (5.14) 
= '^^ {nni;u, (r - n)^) /(r - n) 6{t - r') 

Due to the relation [ToTq)^ — 0,{i — 1,4) only (2,3)- 
components enter the action. The bare values of the 
vertices are equal to 

ArM^,«,^^i;r^) = {^?|^J III (5.15) 
Finally, we write the cubic term in the form 

53 = -2V^ie5i3-y Y1 / dXdXi (5.16) 

ij cr=± 

(t/)^ ^uTalljVx,/?) Bi{X,Xi) (J^adxiTsUjtpXu'y) 

with 

Bnx,X,) (5.17) 
= (imi; ti, u,, (r - n)^) /(r - n) S{t - r'), 

where the bare values of this vertex are 

Br(0,u,«i;rj.) = ||j^j iZl'l (5.18) 

As we will see, in the approximation we consider, a 
generalization for the terms So,Sbi,Sb2 will be necessary 
only in the Id case. 

B. Renormalization scheme 

We use a standard momentum shell renormalization 
group scheme. Separating fast and slow fields in the ac- 
tion we integrate over the fast fields and determine in 
this way the flow of coupling constants as a function of a 
running cutoff. In our case this amounts to a resumma- 
tion of the perturbation theory in the leading logarithmic 
approximation. A quantity y is expanded in a series of 
the form 

y = ^[7ln(...)r a„(7), (5.19) 

n 



and one attempts to find a Taylor expansion of a„(7). 
We assume during the renormalization that the coupling 
constants 7 are small, 7 <C 1. 

In our case it convenient to define fast fields (f> and 
slow fields 'S' with respect to the frequency only. The 
reason is the anisotropy in momentum. As one can see, 
relevant momenta p|| are of the order of lo/vf, while mo- 
menta pj^ do not contribute to the logarithm and enter 
as parameters. Thus we write 

i>{X) = ^{X) + cf,{X), (5.20) 

where the fast fields have the frequencies uj in the in- 
terval, 

KUJc < \^\ < (5.21) 

while the slow ones 'S' carry frequencies 

\uj\ < KUJc, (5.22) 

where Uc is the running cut-off and k < 1. Fast modes 
are integrated over in the Gaussian approximation using 
averages of the form 

(...)o= J d(/>(...)exp(-5o [(/>]). (5.23) 

This results in a change in S 

SS[^] = - In ^cxp ( - 5[* + 0]) - , (5.24) 

that will now be determined explicitly. In diagrams the 
Green's function of the fast modes will be denoted by 
a thick solid line in order to discriminated it from the 
Green's function of slow modes. 



C. Renormalization of interaction amplitudes 

The renormalization of the interaction amplitudes was 
considered in Ref. O Here we merely summarize the 
results, since we will use them later on. Let us note that 
for the renormalization group the symmetric choice a = 
1/2 is most convenient. Relevant diagrams are shown in 
Fig. H 

The result of the analysis in Refilj was that the model 
is reproduced under renormalization and the changes in 
Fi, Ai and Bi can conveniently be written in the form 

5f,; = <Br, (f,;^,;A,),5e (5.25) 

= Q5e, (f,;4;A,)(5C (5.26) 

5A, -<Ba. (f,;6,;A,)(5e (5.27) 

where 

= uuifij^ (^^^ In (5.28) 
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m 




(SA) 



FIG. 4: Relevant diagrams for the renormalization of the in- 
teraction amplitudes as found in Ref. Q 



Here fid in d dimensions is given as 
2 



2, 

4(pFro) ^ 



d = 2 



47r(pi.ro)~2 , d = 3 



(5.29) 



and 



-/(P±) (5.30) 



(27r)'i-i 

Therefore the amplitudes can be written in a scaling form 

r,(0;w,ui;ri) =7,K(0;u,7.i;ri);7O(0)] (5.31) 
Sn^; ^, "i; r^) = K(^; ui; r^); 7°(0)] 

Ai^'''{e-u,u,-v^) = Ar'^Me(e;",«i;ri);7°(^)] 

where 7? = 72 = 7/, 7° = 74 = 7& give the initial 
conditions for the flow. The flow stops at max(6', T/eoo), 
so that 



^{9; u, ui;r±) = ~uuifidf± ( ^ ) In ( max 



T 



(5.32) 



For any perturbative calculation of the spin suscep- 
tibility the amplitudes 7i,/?j^,A^'^ cannot enter. This 
immediately follows from the relation Ait^J-q = —J-q, 
which means that the matrix structure of the Green's 
function in iJ-space becomes trivial for every closed loop 
in perturbation theory. Since in this paper we are in- 
terested in the perturbative sector of the model only we 
give here the relevant RG equations for the backward 
scattering components 



^73 (0 
d^ 



[73(e)]' 



-273 (e)/33+(e) 



dPs (e) 

dC 



(5.33) 



-73(0/33"(0; 



dAg + iO dA+~ (0 



= -2A++ iOlAO-HPtiOy 



d^ 



d( 



-2/33- (0/3,1(0 



(5.34) 
(5.35) 

(5.36) 



There is a subtle point related to the amplitude A3 
Instead of a flow equation the relation 



A3- (073(0 = [/33"(0]' 



(5.37) 



was fixed in Ref. 1 to cancel ultraviolet divergencies that 
would otherwise develop under a change in the cut-off. 
We will come back to this point in Sec. IVII below. 

Appropriate boundary conditions have already been 
specified when introducing the model above. The solu- 
tions of the fiow equations are 



73 (0 = P3 (0 



1 



/33+(0 = A+- (0 = A3-+(0 = 



A++ 
^3 



24* 



where we introduced the notation 

1 



> 0. 



(5.38) 
2,(5.39) 

(5.40) 
(5.41) 



and the backscattering amplitude 7° is defined in 
Eq. 



D. Renormalization of So, Sto, Sti and Sb2 

In this section we consider the renormalization of the 
terms Sq, Sbo, Sbi and 5^2, Eqs. ({0^14.4614.'^ . It will 
be shown that for the one-loop RG considered in this 
paper vertex corrections cancel in dimensions c? = 2, 3 
and, as a result, these terms are not renormalized. This is 
no longer true for d — 1. Unlike in higher dimensions no 
angular integration is performed in one spatial dimension 
and this fact is responsible for the appearing of additional 
logarithmic corrections as will be shown below. When 
selecting the relevant corrections in Id, we have in mind 
a comparison to the well known result of Dzyaloshinskii 
and LarkinS^. 



1. Corrections to So 

This contribution has been noticed before^ but was dis- 
carded, since for the renormalization of the interaction 
amplitudes this term was already beyond the desired ac- 
curacy. The relevant diagram is shown in Fig. [S] In the 
presence of external vertices this term should be consid- 
ered. 
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5S, 



(1) 



ss, 



(2) 



5S, 



(3) 



FIG. 5: Logarithmic correction to So in d = 1. 



After expanding field '^{xi) around a point x one ob- 
tains an expression containing the following integral 



tj: dp 



vpnip 



1 



A{p) (5.42) 



where A is a product of the vertex parts B and cut-off 
functions / appearing in the expression. The crucial 
point is that there is a free integration over the vector 
n.2 and unlike the contributions from diagrams shown in 
Fig. 21 a logarithm can be obtained only in d = 1 . The 
result can be written as 



where 



STZ = u J duiS£,{u,ui) 

X [(1 + n3)/3+/3r + (1 - n3)/33+/33- 



(5.43) 



(5.44) 



It seems natural to interpret this term as a correction to 
Sq. For our purposes it is more convenient, however, not 
to allow Sq to change. This can be achieved by rescaling 
fields ^p after each renormalization step in such a way STZ 
is removed from Sq. This is why we do not write addi- 
tional RG equations here. In turn, this rescaling of the 
fields can lead to additional corrections in the flow equa- 
tions for the interaction amplitudes or external vertices. 
For the interaction amplitudes, it is in fact easily seen 
that taking these corrections into account would be an 
overstepping of accuracy. This is no longer true for the 
external field vertices as will be discussed below. 



2. Corrections to Sb2 

Relevant contributions to the term Sb2 are represented 
in Fig. The correction dsjj^^ is determined by the ver- 
tices A and F and can be written in the form 

(2) I 

SSI2 = 81^7] SaiSj / dXdXi uuiba{x) (5.45) 

X73A+" 6^ ^X,f3^t^^X^.,^ fix - Xi) , 



FIG. 6: Diagrams for the corrections to Sb2- Logarithmic 



corrections from SSj^l^ 



and 5S, 



(2) 



cancel each other in any 
dimension. gives a logarithmic correction in d = 1, but 

for the susceptibility it gives corrections beyond our accuracy. 



where 



and 



++ 



T+KT+ 



f{x - xi) = /(r - ri)(5(T - n) 



(5.46) 



(5.47) 



The correction is logarithmic in any dimension. However, 
the form of is different from that of Sb2 because it 
contains integration over both n and ni, which contrast 
the bare form Sb2, Eq- (|4.48|) Moreover, the matrix K 
breaks the symmetry in g-space (superspace) and, at first 
glance, one should introduce additional renormalization 
coupling constants. 

Fortunately there is another diagram that exactly can- 
cels the previous one. It is also shown in Fig. and its 
contribution equals 



SS, 



(1) 



62 



—8i'r]eaf3-y J dXdXi uuiba{x) (5.' 

^PtP3 5^ i^x.,0^t^^x,n fi^ - ^1) 
In fact, one comes to the exact cancellation 



SS, 



(2) _ 
b2 ~ 



-SS, 



(1) 
62 



by virtue of the relation 

73Ar 



(5.49) 



(5.50) 



that follows immediately from Eqs. (|5.38I5.40|I . 

Finally there is an additional logarithmic contribution 



in Id, (55^2^ , represented in Fig. El It has a similar form as 
SSq, but taking this correction into account would mean 
overstepping the accuracy for our problem. The reason 
is that due to the supersymmetry no diagram for the sus- 
ceptibility can be formed with the help of the vertex Sb2 
without including additional interaction amplitudes. As 
a consequence, the leading correction to the spin suscep- 
tibility in 0? = 1 resulting from this contribution would 
be Sx 7"^ ln(. . .), which is beyond our accuracy. For the 
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A 




FIG. 7: Logarithmic correction to Sto in d = 1. 

same reason the rescaling of the fields, which is necessary 
to bring Sq to its bare form after the renormahzation, 
need not be considered here. This result relies on the 
symmetry in the g-space and this is why it was impor- 
tant to check the cancellation of the terms violating the 
supersynimetry. 

3. Corrections to Sbo 

The analytic expression corresponding to the diagram 
shown in Fig.[7|is quadratic in b but does not contain any 
slow field Therefore, we attribute the corresponding 
contribution to the renormahzation of Sbo, Eq. i|4.46|l . 
Clearly, in this case the rescaling of the fields is not im- 
portant. Since there is a free integration over both the 
vectors ni and n2 for d > I, the correction is logarithmic 
only in d = 1. This is similar to what happens when 
calculating the correction to Sq. 

In d — 1, we write Sba in the form 

Sbo — —-^^11 J dxduidu2h^{x) a{(,,ui,U2) (5.51) 
and set 

cr(e - 0,ui,U2) = 1. (5.52) 
Then the correction 6a to this quantity takes the form 

(5.53) 

4- Correction to Sbi 

There are two separate contributions to Sbi, Eq. 
H4.47(l . represented in Fig. |H1 In d > 1 the slow field 
\E' fixes the vector n in one of the Green's functions G 
only, while the other vector n' is integrated over. As 
a consequence, a logarithmic correction is obtained only 
in d = 1. In one dimension the rescaling of the fields 
(cf. Sec. IV DT)) is also important and gives an additional 
contribution. 




SSbi : 5fj,a + Sfib SSbi: Sfic 



FIG. 8: Logarithmic corrections to SSbi in d = 1. 



Considering the case d = 1 we present Sbi in the form 
iSbi = \J —2ivr\ J dxdzidz2 ba{x) 

'X-^0T3l^{^,Ui,U2)lpa{x,U2) . (5.54) 

Trivial " angular" integration in Id (weighted summation 
over the directions) has been performed. The operator 
V is defined as 



V{^,Ui,U2) = 




Here /i^ = /ij(^, ui, U2), and initially /ii(^ — 0,ui,U2) = 
1, i — a,b, c. In this case 

V{^ = 0,ui,u2) = d{x,u2) (5.56) 

and we come back to the original form displayed in 
Eq. PTzjl . 

The diagrams in Fig.|Hlrepresent corrections to fii. The 
left diagram determines corrections 5fia and Sfib 

6fia = -y (/33" A+- + /33+A3— ) 5^ - -5fib (5.57) 

The contribution S^ic consists of two parts 

= (5.58) 

The correction Sfii^^ is represented by the right diagram 
in Fig.lHland reads 

S^ii'^ - ~U2 {P+A-+ + P^A+-) SC (5.59) 

The correction Sfi^^ is due to the rescaling of the fields 
has to be performed at each RG step to keep the form of 
Sq fixed 

Sn<;?^ = -2u2l3+l3^S^ (5.60) 

Note that the forward scattering components drop out as 
could be expected. 
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5. RG equations and their solution 



We found logarithmic corrections to the vertices Sbo^ 
Sbi and Si,2 in dimensionahty d = 1 only. This means 
that these terms are not renormalized in the first order 
in the dimensionalities d = 2, 3 and the vertices 7, (3 and 
A given by Eqs. (|5.38I5.4U|I are sufficient to determine 
the susceptibitlity. 

At the same time, the renormalization of the vertices 
SbOi ^bi and 5^2 is very important in d = 1. Both func- 
tions fii and cr from Eqs. 1)5.511 I5.55|l do not have a 
simple form and one should write and solve proper RG 
equations. For the function a related to 6Sbo, we write 
a = (7(^,1*1,^2) and using the correction da, Eq. H5.53)l . 
obtain the following differential equation 



da 



1 f 6^* 2^* 



(5.61) 



With the boundary condition a{^ = 0, ^1,^2) = 1 we 
obtain 



0'(C,Ml,U2) 



(5.62) 



1 + -U1U2 Tj 



2C 



The corresponding differential equations for /^^ are to be 
obtained from the forms of the corrections, Eqs. 1)5. 57t 
15.601). and can be written as 



-Ul 



-U2 



*\3 



i^ + Cb) 

34* 



(5.63) 



(5.64) 



with the boundary conditions /ii(^ = 0) = 1. Integrating 
these equations we obtain (only fia and fic will enter our 
results) 



Ma 



Cb 



2 l,(e + C)' 
,2 



1 



(5.65) 



24* 
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The calculations presented in this subsection allowed 
us to obtain all effective vertices entering the RG scheme. 
This gives us a possibility to calculate the susceptibility 
for all dimensions d = 1,2,3. The result for d = 1 is 
well known"^^ from a renormalization group treatment for 
the initial electron model. We will reproduce now this 
result using the derived equations in order to check the 
formalism of the bosonization used here. Only after that 
we will concentrate on calculating the susceptibility in 
the higher dimensionalities d = 2, 3. 



VI. SPIN SUSCEPTIBILITY IN d = 1 

We can now determine the temperature dependent cor- 
rection to the static spin susceptibility in d = 1. 



SX2 



A 



FIG. 9: These two diagrams with renormalized external ver- 
tices (5x2) and interaction amplitude (^Xs) determine the 
correction to the susceptibility in d = 1 together with the 
diagram of Fig. Q 



We write the correction to the susceptibility as 
Sx = <5xi + Sx2 + Sx3 



(6.1) 



and denote the contribution from SSbo, Fig. as Sxi- 
The second term i5x2 in Eq. I|6.1|l is the contribution 
that corresponds to the diagram shown in Fig. |51 on the 
left hand side. 

Here only the renormalized vertex of Sbi enters but no 
additional interaction amplitude. A diagram with this 
property does not exist for vertex Sb2- Finally, the di- 
agram shown on the right hand side of Fig. |3 gives a 
correction termed 5x3 ■ It involves a renormalized inter- 
action amplitude. The corresponding expressions take 
the form 



Sxi = 



J duidu2 a{^,ui,U2) 



(6.2) 



6X2 ^ 2j^y duidU2dU3 U2 ^J.ci^,U3,U2)^J.a{^,Ul,U2) 

(6.3) 

6x3 ^ J duidu2UiU2 ^ A3'^(^,iti,-U2) (6.4) 

Q,/3 = ± 

We do not write factors of 77 in d = 1, since we want 
avoid unnecessary complications while focusing on the 
leading temperature dependent corrections. In Eq. 1)6. 4|l 
we introduced the interaction amplitudes Ag^. Naively 
one would expect amplitudes Ag'^ Eqs. H5.38I5.40|) to 
enter here but this would not be correct. In fact, this 
question is intimately related a subtle point related to the 
renormalization of Ag"^ already alluded to in Sec. IV Gl 

When calculating corrections to A;^~ within the renor- 
malization scheme, the authors of Ref 1 found ultraviolet 
divergencies that could be cancelled only provided the 
condition A;^~73 — {(3^) is imposed. Since and 
73 can be determined independently, this condition fixes 
A;^~ . For large temperatures, where one can use the bare 
values of these amplitudes, this relation is automatically 
fulfilled. It is crucial to note now that it was necessary to 
fix A;^~ only because this amplitude itself enters ^2 [</>], 
where fields (p are the fast modes. 



17 



Returning to the diagram for (5x3, Fig.|51we see that the 
frequencies flowing through the Green's functions are de- 
termined by the external vertices and therefore are van- 
ishingly small. In particular, they are smaller than any 
frequency considered in the renormalization scheme. We 
argue that the part of ^2 that contains the fields at van- 
ishingly small frequencies should be split off from the be- 
ginning and when separating fast and slow modes, it must 
always contain slow fields only. Correspondingly, the in- 
teraction amplitudes, termed Ag'^, are renormalized but 
do not play any role when calculating corrections to the 
interaction vertices. In such a situation, there is no rea- 
son to fix as was done previously for Ag ~. Instead, 
one should follow the renormalization group scheme and 
derive a proper RG equation for A;^~. 

The relevant diagrams have been already presented in 
Fig^and the result of the RG procedure can be expressed 
by the equation 



dA 



(6.5) 



The solution of Eq. H6.5|l with the initial condition 
A3""(^ = 0) = takes the form 

2 1 



A- 



(6.6) 



This should be contrasted with 



A,— = 



1 



(6.7) 



We checked our reasoning by a perturbative calculation 
at order 7^ , where the difference between A3 
is already noticeable. 

Finally, we use the identities 



ig and A3 



1 + X 



duidu2 U1U2 {zi2 + + 22^2) (6.8) 

= f duidu2 u\ [zl^ + 22^2) (6-9) 

where Z12 — 1/(1-1- U1U2X), that can be checked by a 
direct computation of the integrals. Then, recalling that 
X — 2^bh\{£F /T) we come to the following temperature 
dependent correction to the spin susceptibility. 



5x{T) = 



l + 27bln(eF/T) 



(6.10) 



This result has first been obtained by Dzyaloshinskii and 
Larkin'^*. Eq. H6.1()|l serves as a good check of the 
bosonization approach used here. Actually, the calcu- 
lations within the framework of the bosonization method 
of Ref.^ are most difficult in d = 1. It is clear that 
this method is less convenient for calculations in Id than 
the other well developed ones^S. However, calculations 
in d = 2,3 are somewhat less involved and the present 
approach is the most convenient tool for calculations in 
these dimensionalities. In the next Section we concen- 
trate on such calculations. 



VII. NON-ANALYTIC CORRECTIONS TO 
SPIN SUSCEPTIBILITY IN d = 2, 3 

Non-analytic corrections to the spin suscep- 
tibility have been considered in several works 
befor3^^'^i^''^'^'*i^''^'''^i^7ii^i22. A linear in T behavior 
at order 7^ was obtained in 2c?, while the potential 
analog in 3d, a InT behavior, was found to be absent 
and the first correction in 3d was proportional to T^. 
We will show now that there are logarithmic corrections 
to these results and sum up the leading logarithms. 

Let us repeat that, as it has been demonstrated in 
Sec. lVDI the terms Sba, Sbi and Sb2, Eqs. (|4.46I4.48|I . are 
not renormalized in dimensions d > 1 . Therefore we can 
perform a perturbative analysis with the renormalized 
interaction vertices 7, /9, and A, Eqs. H5.38I5.40|) . keeping 
the bare values of Sbo, Sbi and 5^2. 

The relevant diagrams leading to T'^~^ corrections are 
displayed in Fig llOl The solid lines carry the frequencies 
Lo and the momenta k of the order of T and T /vp, re- 
spectively. They are smaller than characteristic energies 
in the Green function entering the vertices because the 
latter are responsible for the logarithmic contributions. 
This means that the vertices can be taken at zero ex- 
ternal frequencies and momenta and this is the reason 
why one may just take the values of the vertices from 
Eqs. (|5.38I5.40|I . The same procedure has been used in 
Refi for calculation of the specific heat. Putting the bare 
values for the vertices 7, j3 and A would give the petur- 
bative results of Refs.ii*i2iiii4 in d = 2, 3. In this limit, 
the diagrams of Fig. ^] correspond to the conventional 
diagrams considered in those works. 

As concerns diagrams containing the amplitude 7/ of 
the forward scattering, we did not find any logarithmic 
contributions. This is because one obtains integrals of 
products of Green functions containing poles in the same 
half plane of complex variables kn. 

Considering the contributions of the diagrams in Fig. 
1101 and comparing them with self-energy and vertex cor- 
rections in Id, Fig. Elone can see that there is a close anal- 
ogy between the terms responsible for the logarithmic 
corrections in one dimension and those responsible for 
the non-analytic behavior in higher dimensions. Within 
our formalism, the main difference between the two cases 
is the additional angular integration in dimensions d > 1. 

Let us now turn to the computation of the diagrams in 
two and three dimensions. Calculating the terms of the 
perturbation theory corresponding to the diagrams dis- 
played in Fig. 1101 one finds that some of them show un- 
physical divergencies in the limit of vanishing momenta 
and frequencies. Therefore one should sum up certain 
diagrams first before taking the limit. 

To demonstrate this feature explicitly, let us consider 
the backscattering contribution for the diagram 3 of Fig. 
[TUl When evaluating the term X3(q, £) corresponding to 
this diagram one finds terms containing the product 



^ni(q,£)5n2(q, -e) 



(7.1) 
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are included in Appendix lEl. Using the mutual relations 



(7) 



FIG. 10: These are the diagrams responsible for the non- 
analytic temperature dependence of the spin susceptibility. 



where q and e are the external momentum and frequency. 

The integral over the internal momenta is ultraviolet 
divergent and must be cut with the help of the function /, 
Eq. H3.6|l . At the same time, the limits q — + 0, e ^ in 
the diagram 3, Fig. 1101 cannot be taken unambiguously. 

In order to get rid of such unphysical divergencies, we 
note that this term contains the product of the inter- 
action amplitudes Aj^~r3. A closer inspection reveals 
that this term is intimately related to the renormaliza- 
tion of A and the ultraviolet divergence encountered dur- 
ing the renormalization (see Fig. 0)). It follows from the 
results of the renormalization that this divergence must 
be cancelled with the help of diagram 4 using the relation 

A3~73 = (/33")'- 

The relation to the conventional perturbation theory, 
which is made obvious by including the dotted lines into 
the diagrams, in fact strongly suggests first to group sev- 
eral different diagrams before evaluating them. These are 
the groups 



Xa = Xl + X3 + X4 + X5 
Xb = X2 + X6 + X7 



(7.2) 
(7.3) 



Strictly speaking, diagram 6 differs topologically from 
diagrams 2 and 7. It nevertheless turns out to be ad- 
vantageous to combine them, since the expressions are 
similar at low energies that are considered in the model 
we use. 

For convenience of the reader, and since the algebra 
is rather tedious, explicit expressions for the diagrams 



between the seven interaction amplitudes A^'^, and 
73, one finds rather simple expressions for Xa and Xb, 
that allow to take easily the limit q ^0, e ^ 0. After 
taking this limit they read 



Xa 



-327?2r^ J dp J dnidna (7.4) 



Xb 

where 

Y{p,( 
and 



{vpT^iP — iw)^{vFn2P + i^Y 
2,2TfT^ i dp j dnidn2 

iu!{vF'n.2P " i^) 



F(p,fMn5) 



(7.5) 



(u_Fnip - iu!)^{vFn2P + iuj] 

ib{e) 



d'^re-'P'"- 



1 



F(p, nin2) 



/(r) 



(7.6) 



X{e) = ^^ldlb{e) in (max{(?, r/£o}) (7.7) 

The numerical coefficient fid was introduced in Eq. H5.29|l . 
The integration over ui, U2 was performed with the help 
of the following relations 

j duidu2u\u\(^ziZ2 + 2{z\z2 + Z1Z2) (7.8) 

+ {ZIZ2 + Z^zl) + Izlzi) = \ 

J (1 + a::i)(l + a:;2) 
j duidu2u\u2 (2zfz2 + 2{zfz2 + zizl) (7.9) 

+(^'^^ + ^^^^')) = (i + .0(i + x2) 

where 

Zi = — i (7.10) 

1 + UiU2Xi 

The non-analytic contribution to the spin susceptibil- 
ity is found from the small region of phase space, for 
which the angles rii and n2 are close to each other, 
|ni — ^ 1. We therefore introduce 

n = (ni -h n2)/2, m = rii - n2 (7.11) 
P\\ = Pn, p± = p-_p||n (7.12) 

and perform the integration in p|| in Eqs. (|7.4() and (|7.5|l . 

As a result, we obtain the following formula for the 
non-analytic correction to the spin-susceptibility in di- 
mension d > I 



6x{T) = Xa{T) + Xb{T) 



(7.13) 



256 2^Y-| 13 f d^-W f, , 



3(wi?mp^)^ — iuj 



77 X2 I A 2^3 ^(^'ll ~ P-L' 
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The main contribution to the integrals in Eq. (|7.41 17. 5|) 
comes from p|| of the order of T/vp and this is why we 
can set p|| ^ in the argument of the function Y . 

Equation 17 . 131 contains a sum over bosonic Matsubara 
frequencies and we write this sum symbohcahy as 



(7.14) 



Technically it is more convenient to calculate the devia- 
tion from the zero-temperature limit instead of comput- 
ing the sum, i.e. to calculate the quantity 



5x{T) = 5x{T) - 6x{T = 0) 



(7.15) 



Using the Poisson formula the temperature dependent 
correction to the susceptibility 6x{T) can be represented 



. Ilu 



(7.16) 
(7.17) 



The further evaluation is slightly different in dimen- 
sions d = 2 and d — i and we discuss the two cases 
separately. 



A. Non-analytic correction in two dimensions 

Rescaling the momentum and integrating over the an- 
gle n we obtain 



5x{T) 




(7.18) 



m 



The integral over |m| is logarithmic and therefore not 
very sensitive to the upper limit that can safely be set 
to 1. One notices that the momentum dependence of Y 
is crucial here coupling the integrals in |m| and k. liY 
were independent of the momentum, the k integral would 
be equal to zero, whereas, at the same time, the integral 
over m would diverge at the lower limit. 

Fortunately, this uncertainty can easily be avoided tak- 
ing into the momentum dependence of the function Y . 

After introducing the Fourier transform of Y 



Y{\rVA 



2tt 



Yip^,p»^0,d) (7.19) 



the momentum integration can be performed with the 
help of the identity 



dk 3fc2 



1 



ikb 



(27r) (1-Hfc2)3 



2„-|6| 



(7.20) 



where b = 2r|u;|/(wi?|m|). 

Then, we use Eq. (|7.17(l and obtain the following ex- 
pression 




2,,,2„-/3|c^| _ 



= T\2x 



, cothx 
sinh^ X 



(7.21) 



where x = 27rT|r|/(t;i?|m|) = ttT/?, and change the inte- 
gration variable from |m| to x. As a result, we find 



SxiT) = - 



AT 



drl2{a)Y 



where 



Ma) = 



dx 



/ 2x coth X 2 
\ sinh^ X a;2 



2nT r 1 

£oo To X , 

(7.22) 



(7.23) 



and 



2-nT r 1 
Eoo ro A(f> ' 

We reintroduced formally an upper cut-off for the 
integration over |m| but it will drop out from the final 
result. One can see that the essential r as controlled by 
the function Y (and thus /) are small, r < tq, while 
essential x in the integral are large, x > 1. This means 
that the main contribution in the integral over the angles 
6 comes from 9 of the order of T/soo- Therefore we can 
with logarithmic accuracy set 9 — in the argument of 
Y. In turn, it means that the integral over x is rather 
insensitive to the lower bound as long as a; < 1 and we 
may safely extend the integration range in x to the in- 
terval (0, oo). Then, the integrations over x and r can be 
easily performed. To this end we note that l2(0) = — 1 
and introduce the notation 



Y{9) ^Y{9, p = 0) 
to formulate our result for the susceptibility dx, 

5x^\T)^2^^—xl''Y{9^Q) 



(7.24) 



(7.25) 



In Eq. ((7123, xl"^ = rnj-K. The vertex part ?], Eq. Mj^ . 
should be understood as T]a=\ and it is a result of an 
additional summation of ladder diagrams including 52, 
in close analogy to the discussion in Sec. lIVFl The limit 
= corresponds to the backward scattering. Before 
further discussing this result in Sec. lVIl'Cl we turn to the 
three dimensional case. 



B. Non-analytic correction in three dimensions 

In 3(i one obtains from Eq. (|7.13|) after rescaling of mo- 
menta and integration over |n| the following expression 



W) 



327r 



T 




1 d|m| 1 
2-K ImJ 
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[ dPk 3(emk)2 - 1 / 2|w|k , a 

J (27r)^ ((emk)^ + l)'^ V VF\m\ 7 

(7.26) 

where em = m/|m|. 

We see that the integral over m in Eq. H7.26|l is loga- 
rithmic. However, corrections of the form 5x oc 7^r^ InT 
are absent and this is due to the fact that the integral 
over k vanishes provided the momentum dependence of 
Y is neglected. Nevertheless, if the function Y depends 
on the momentum p_L the entire integral is finite. 

After introducing the Fourier transform of Y as 

Y{\vie)= I ^e^P-ry(p||^o,px,0) (7.27) 

it is convenient to decompose the vectors k, p± , r into 
components parallel and perpendicular to Cm, such that 
r = (f|| ,r±). Then, one can then proceed in close analogy 
to the calculation in 2d. Using Eq. (|7!^ and Eq. (|7!^ 
one arrives at 

6x{T) = -v^ / d\Ma)Y( rl0= U_ , 

Vp J V (Too ro xl 

(7.28) 

where 

^ , . r , / cothx 2 \ 

and a = ^^-[^s^- Again, we reintroduced the upper 
cut-off A(/) for the m integration. 

The integral over x in Eq. H7.28|l shows a some- 
what stronger dependence on the lower cutoff a as com- 
pared to the two-dimensional case with 23(0) = —1/3, 
2^3(1) ~ —0.28. This means that angles larger than T/soo 
start contributing more significantly. Still, the dominant 
contribution to the integral comes from x > 1, so that 
we can set 6* = in the argument of Y with logarithmic 
accuracy. Then, we come to the following result for the 
temperature dependent correction to the spin suscepti- 
bility 

5x'\T)^^Tf'^xl''Y{e^Q) (7.30) 

o £ p 

where Xo'' = tnypj-K'^ and Y(&) ~ y(p = 0,6*) and r\ — 
T?Q=i, see the remarks below Eq. (|7.25|) . 

We see that, as in the 2d case, this correction is deter- 
mined completely by the backward scattering (6* = 0) . 

C. Final results in d = 2, 3 

1. General results 

The quantity Y{9 — 0) can be considered as the square 
of an effective temperature dependent backward scat- 
tering amplitude jb{T), and we write it in the form 



Y{0 = 0) =72(r), where 

f d^-^v^ 7-L (ft) 

lb{T)=ib —j-i _ / / , (7-31) 

^0 1 + /^(£J^)X(T) 

76 = 76 (^ = 0), X{T) - /irf76ln(eoo/r) and 7±{r±/ro) 
was defined in Eq. H5.30|l . 

So, our results can be written in the most general form 

as 

Sx^'iT) = 2rj'jUT)—xl' (7.32) 
Sx^'iT) = f^rj'j^^{T)^ xl" (7.33) 

and we remind the reader that 77 = 770=1, where 77 is 
determined by Eq. H4.49|l . 

If we replaced jb{T) in Eq. (|7.32() by the bare cou- 
pling constant 76 for d — 2, we would obtain the previ- 
ously reported linear T dependence of the non-analytic 
correctionsiiiiSiiii^. This replacement means neglecting 
the rcnormalization of the interaction constants discussed 
in Sec IV CI If we set the function 7b (T) equal to the bare 
value 7fc in d = 3 we would obtain the correction 6x'^'^{T) 
proportional to T^, which is regular in T^. This means 
that the first non- analytical T^ln(eoo/T) term in 3d is 
of the order 7^ . 

In the limit of small X {T) ^ 1 the temperature de- 
pendence of 7j (r) takes the form 

lUT)-lb~2jlc,ln'-^, X(r)«l, (7.34) 

where 

c. = M./^7l(^) (7.35) 

The factor Cd depends on the precise form of the cut-off 
and can be estimated only. It is roughly of the order of 
unity. Eq. H7.34|l shows that the first logarithmic in tem- 
perature corrections contain the prefactor 7^^ both in 2 
and 3 dimensions. This rather high order in the coupling 
constant 7;, is, apparently, the reason why the logarith- 
mic corrections to the susceptibility have not been no- 
ticed previously in the diagrammatic expansionsii^i^iiii^ 
(see, however;^ for 2d). 

In the limit of large X (T) ^ 1 one finds the following 
asymptotic temperature dependence of 7^ (T) 

7,^(T)(x(ln^)"', X(T)»1. (7.36) 

More explicit formulae can only be written using a model 
cut-off function and this will be done in the next section. 
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2. Results for a model cut-off function 

We choose the fohowing model cut-off function 
A (-i/ro): 



?0 



— ^ exp I 



(7.37) 



where f2d_i is the d — 1 dimensional solid angle. 

Performing the remaining integration for this case one 
obtains the following temperature dependence for the ef- 
fective backward scattering constants 



27bln[l + X(r)/2] 
X{T) 

2^7bL^2 [-X{T)/2t:] 
X{T) 



(7.38) 
(7.39) 



where lji2{x) = Tl'^^i^^ 1^"^ the polylogarithm func- 
tion. 

In the limit of small X (T) <C 1 the temperature de- 
pendence of the susceptibility computed with the model 
cut-off function takes the form 



Sx^" (T) 



2^^,^^X^'^(l-27.1n^) (7.40) 
^^(l->ln^) (7.41) 



_r,2 2::_ 3d 

O bp 



where we put r'^^ ^ pp for simplicity. It should be 
stressed once again that the coefficient of the logarith- 
mic correction cannot be determined rigorously within 
our model. 

In the opposite limit of very low temperatures, 
X (T) ^ 1, asymptotic expressions for the corrections 
to the susceptibility can be written using the model cut- 
off function of Eq. (fTTzjl 



1 2T.2dln' (476 In ^) 



2^ 



In 



2 e„ 



^x'^iT) = -V'—Xl' 



3^1n4(47r7bln^ 



In 



2 



(7.42) 
(7.43) 



Again we used Tq ^ ^ Pf for simplicity. The asymptotic 
behavior 1/ (ln(eoo/r))^ in these equations is not very 
sensitive to the form of the function / (k) as can be seen 
from Eq. lf73T|l . 



VIII. DISCUSSION 

We have calculated non-analytical logarithmic in tem- 
perature contributions to the spin susceptibility of a d- 
dimensional electron gas for d — 1,2,3. We used the 
bosonization method recently developed in Refi and 
demonstrated that it can give results not only for the 
specific heat as in Ref.^ but also for the spin susceptibil- 
ity. 



The main contribution to this quantity comes from ef- 
fective spin modes that interact with each other, which 
leads to the non-analytic logarithmic contributions. Al- 
though we consider isotropic systems, the low tempera- 
ture behavior is determined by spin excitations moving 
antiparallel to each other. As a result, the non-analytic 
contributions are determined by the backward scattering 
showing that there are one dimensional features also in 
the dimensions d = 2, 3. 

The final form of the temperature corrections to the 
susceptibility in two and three dimensions is given by 
Eqs. (|7. 3117. 33(1 . Although in 2d the correction to the 
susceptibility x is very similar to the correction to the 
quantity C (T) /T, where C (T) is the specific heati, they 
are quite different in Sd. The first logarithmic contribu- 
tion to C (T) /T is of the o rder 7^, which is a well known 
result for 3d (Refs. Ulsll^ . At the same time, the 
expansion of the susceptibility in the logarithms starts 
with the term of the order of 7;^, which shows that the 
non-analytical in temperature corrections exist for this 
quantity in three dimensions, too. 

Using the bosonization scheme of Refi we have also 
reproduced the temperature dependent correction in one 
dimension, Eq. (|6.1U|) , that has been obtained long ago^ 
using a renormalization group approach for the initial 
electron model. 

The temperature dependent correction to the suscep- 
tibility in 2d was calculated recently by Shekhter and 
Finkelstein^ using direct diagrammatic expansions for 
the initial electron model. In the approach of Ref4£, 
which was tailored for the calculation of the spin sus- 
ceptibility in d = 2, the renormalization of the effec- 
tive backward scattering amplitude is attributed to all 
Cooper channel harmonics, while no cut-off function was 
used. In the formalism of Ref. 1 which we studied no 
decoupling in the Cooper channel is introduced in addi- 
tion to the particle-hole channel in order to avoid over- 
counting in the region of phase space close to backward 
scattering, which turned out to be most important (For a 
more detailed discussion of the role of the Cooper chan- 
nel in the bosonization approach see Sec. VII C of Ref. 
1). In fact, the renormalization of the backward scatter- 
ing amplitude is obtained in this way as well, non-zero 
angular harmonics are however not included. 

It is important to mention that in some cases not all 
non-analytical corrections are accounted for by the back- 
ward scattering. Interesting contributions of the type 
T^lnT to, e.g. specific heat in three dimensions, are 
given by three-loop diagrams in the language of the elec- 
tronic Green functions (Ref. Jj) and they cannot be ex- 
pressed in terms of the backward or forward scattering. 
However, these corrections are proportional to higher 
powers of the interaction constant and are smaller than 
those given by the backward scattering unless the tem- 
perature is very low. In the latter regime the effective 
backscattering amplitude is very small, Eq. H7.36|l . and 
the contribution of three and more loop diagrams can 
become the most important one. Contributions that are 
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not reduced to the backward scattering was discussed in 
Ref. 113 for the spin susceptibility in two dimensions. 



APPENDIX A: DERIVATION OF FORMULA 



EQ. l|4:29jl FOR Zs[h] BY EXPLICIT 
CONSTRUCTION 



Acknowledgments 

We acknowledge stimulating discussions with I.L. 
Aleiner. We thank A. Shekhter and A. M. Finkel'stein for 
bringing to our attention their results^ prior to publica- 
tion. G.S. is thankful to the Feinberg Graduate School, 
to the grant Transregio 12 "Symmetries and Universality 
in Mesoscopic Systems" and GRK 384 of the DFG for 
financial support. The work was completed during the 
visit of K.B.E. to the Weizmann Institute who acknowl- 
edges a support of the EU Transnational Access program 
and the hospitality of the Institute. 



In this appendix we explicitly construct the supersym- 
metric representation of Zs [h] in Eq. (|4.29|) . 

It is a straightforward application of the results of 
Ref. that [h] of Eq. H4.1() in the main text can be 
rewritten as 



Z,[h] = exp (Si^^ uhnix) L-x\, dr'Ki^')^ 



(Al) 



where the form of L^^, will be specified in the following. 
Here X = (r, t, n, u) and the hat in indicates that 
integration in this formula is over the full solid angle. 



a,f3 



Here a, f3 are spin indices and averaging is defined as 

(...)=y'l?(v,cp)(...)e-^[^''^l (A3) 

where 



X 



M 



Vx 



(A4) 



and supervector tp has been defined in Eqs. (|4.14|) . 14.15|l . 



ip = v^^A, A = cr. 



(A5) 



The fermionic part of tp, ip takes care of the normaliza- 
tion via identity 



I?(S,S)e-^[«'«l 



Finally 



AMn 



-iL'l -L' 



(A6) 



(A7) 



TCI / H 



We repeat that L' = {L + L^)/2 and L" = ~i{L-D)/2 
are hermitian. The explicit form of KM is 



AA/„ — ivonW — Aidr + 2iuhr, 



(A8) 



The matrix Ai acts in H space and is written in 
Eq. H4.24|l . Restricting the angular integration to just 
one half sphere we can cast formula Eq. (|A4(I in a new 
form by introducing supervector <p as 



ip{-n) 
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C has to be modified accordingly. 



where 



M„ 

A/_„ 



(/.^ (AlO) 



(All) 



The expUcit form of AM is 

AMn = iDonSsV - kidr + 2mH„ 



(A12) 



Here we introduced S3 = 0-3"^ and Hn(a:;) of Ea._ (|4.25l) . 

Finally the number of field components in ^, 4> is dou- 
bled once more by introducing the electron-hole {eh) sec- 
tor. This can be done by introducing new vector 



1 / 0* 



tp = ip^A (A13) 



Now 



Jx 



M + iS 



where 



Mn 



Mn 

M^' 



cjyx (A14) 



(A15) 



TR 



Note that the transposition for M"'" includes derivatives. 
Here matrix T3 = cr^"'^^ acts in eh space. The explicit 
form is 

AMn,u = -iwoTaSgnV - Aidr - 2iT3]H„ (A16) 

Now we can write an appropriate generalization of 
Eq. (|A1|I . Eq. (|A3p and make contact to formula 
Eq. H4.29|l in the main text. We write 



C[ip,ip]=-i2iy xPxi^ + iS^)'>Px (A17) 
Jx 

where Ti, — AAi (factor of 2i/ is introduced for conve- 
nience) and the averaging with respect to this Lagrangian 
is defined as 

(...) ^ y V{xP, ^) (. . .) e-^I'^'^l (A18) 

Using 



^h(^) = dxM^ix)J^o, ^h(X) = {CJ'hixyy (A19) 
where !Fq is defined in Eq. (|4.25|) and 



dxia) 



u[adr + (1 - a)wonVS3] ) ^^'^^^ 



one verifies that 



(A21) 



(an(a)h„ {Si* - Si*) + h„ {Si + Si) + (n ^ -n 
where 

dn{a) = u{adr + (1 - a)ivonV) (A22) 
Using l^S]^ x^p*x''^ <5n,n"^u,u' one obtains 

(x) 



I[h] = [ T^{X) (i^xipx') W) = \ I hi 

J XX' ^ J XX' 

{{S^x + Six) {Sl:x - Stx)) On{a)hi{x') 

(A23) 

Summation over spin indices 7, /3 is implied. Using fur- 
ther Eq. (jA2|) one finds 



X[h] = ^ J dudx{dn) h^ix) (uL^\^ (drhn) (x)^ 



(A24) 

The last line can be simplified by noticing 

L-iJwonV - dr)h^ix) = h„(x) (A25) 

This equality holds, since ft,h = h x h = 0. The result is 
Zs[h] exp (^-4^1.2 y ^ KiX) (iPx^) Fi,{X')^ 

X exp -oi) (2;)^ (A26) 

This formula is used in the main text, Eq. 

APPENDIX B: DIAGRAMS OF FIClTHl 

In this appendix we give analytic expressions for the 
diagrams displayed in Fig. ^| We introduce 

A^^^^(zi,Z2,k) (Bl) 
^dj.g-*r^.i.. (SliTa, wi, U2, ri)7(r) 

for A and use similar notation for amplitudes B and F. 
To simplify expressions let us write 

4 

A^(zi,Z2,k)=5^ J2 Ap''^{z,,Z2,k){j^,,T7S 

J— 1 <7k=± 

(B2) 
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and suppress the trivial dependence on (zi, 22)- Then the 
relevant expressions read 

Sxi{q,£) = -4r/2T^ / dpdzidz2 u\u\ (B3) 

xstr(A^(p)TW (p, g)A^(p - q)Ti'^^ (p, 9)^ 



^X2(q,e) = -Sr/^T^ / dpdzidz2 u?M2 (B4) 
xstr( A^(p)t(2'') (p, g) A^(p)T(f (p)) 



4 /. 

<5x3(q,e) = 8772^ dpdzidza "?M2 (B5) 



'i— 1 a; 



;r,(p - q) J'oT^^Hp^ 9)A>(p)Ti'^(p, g)^o 



4 /■ 

^X4(q,e) = / dpdzidz2 ufi 



i—l u 



xSr(p-q)Sr(P) (B6) 
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1 "2 



X5(q,e) = -16??'f]^T^ /" dpdzidza 

2 — 1 CTi CJ 

X 6- (P - q) TTlT^^:^ (P, g) A>(p)T(f^'') (p, 

(B7) 



X6(q,e) = -16r/2^^r^ J dpdzidz2 uiu^ 

(71 i — l LO 

Br (P - q) :^TLt' A>(p)T W (p, g)^o (B8) 



X7(q, e) = 8?7^ ^ ^ / dpdzidz2 u\u2 

<Tl,<T2=± i=l 

(B9) 



where 

= [(ic^ + vipE3)-r+(ie + ViqS3)] 

(BIO) 

T(f°)(p,g) = r3(zc^ + VipE3)e^,b)S„ib + g) (Bll) 

b) - T3(zc^ + V2pE3)en.(p) (B12) 

T+{ie + viql]3) + T_(iu; + vipl]3) 
xg„,(p)^„,(<z) (B13) 

(ze + viqi;3) - T_ (icj + V1PE3) 

x^„,(9)g„,(g-p) (B14) 
r+(iw + V2PS3) + T_(ie + V2qE3) 

xe„,b)e„,(q) (B15) 
(iw + vipi;3) - T+{ie + viqSs) 

xg„,(p-g)gn,(p) (B16) 
T+{itjj + V2PE3) + T_(i£ + V2qCT3) 
x^„,(p)g„,(q) (B17) 
= T3(ic^ + vinipE3)eni(p) (B18) 
T+{itjj + V2PE3) + T_(i£ + V2qS3) 
XT3g„,b)^„,b + g)^„,(g) (B19) 
T(J»)(p,g) = r3(ze + vipl]3)e^,(q)e„,(p + <z) (B20) 
Ti^'Hp) = T3(zc. + V2pl]3)e„.b) (B21) 

Four dimensional notation was used p — (w,p), q — 

(£,q)- 
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